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Abstract. 

It is known that the quantized vortices in a supcrfluid can be described by a 
dual electromagnetic model through the duality transformation. Recently a new 
technique, which can selectively remove atoms from a Bose-Einstein condensate, 
was applied to the vortex lattice state of the rapid rotating Boson gases. The 
increase of vortex number, dynamical formation of giant vortex, the oscillation 
and damping of vortex motion were observed in these experiments. In this paper 
we will discuss these observations in the framework of dual description and show 
how to understood these observations naturally. 
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1. Introduction 

It is well known that quantized vortices play an important role in understanding 
many phenomena of superfluid, superconductivity and even astrophysics. In the latest 
several years a remarkable experimental development has been made in the rapidly 
rotating Bose-Einstein condensate(BEC). A large amount of angular momentum is 
imparted into the system and its rotation frequency becomes faster and faster. The 
ground state of the fast rotating BEC was observed to be a regular triangular vortex 
lattice at nearly zero temperature. This provides a new stage of studying both the 
equilibrium and the dynamic properties of the vortex matter. 

Recently two novel experiments were reported by Cornell's group in JILA. 1 2 
Atoms in the hyperfine spin \F = 1) state, were rotated rapidly by the evaporative 
spin-up mechanism. About 180 vortices were observed in the cloud and they were 
arranged into a regular lattice structure. And then, an atom removal laser, whose 
frequency was tuned to the transition energy between the \F = 1) state and the 
\F = 0) state, was applied to the center of the condensate, and the recoil from a 
spontaneously scatted photon blasts atoms out of the condensate. In this way, the 
atoms with lower angular momentum were selectively removed from the system. 

Many exotic features were observed when this technique is applied. After the 
atom removal laser was applied continuously, the vortex number contained in the 
cloud increased form 180 to about 230, and the vortices began to aggregate to the 
center of the condensate, eventually they formed a giant vortex. If the laser was 
applied for a short limited time, it was reported in Ref.p^ that a periodic formation 
and disappearance of the giant vortex core was observed after applying a strong atom 
removal pulse, and in Ref.[2j that the Tkachenko oscillation of the vortex lattice would 
be found after a weak pulse. Both kinds oscillation were found damped. PP0 

A number of questions are raised from these experiments, such as how to 
understand the vortex number increase and the vortices aggregating to the center 
of condensate, why the strong and weak laser strength will cause two different kinds 
of oscillations, how to understand the relation between the oscillation frequency and 
the initial conditions of the condensate, and which mechanism leads to the damping. 
We notice that recently there are lots of theoretical efforts around these experiments, 
some of them try to understand these observations numerically;!?,, and some of them 
focus on the Tkachenko wave of the vortex lattice stated [5] @ HI- However, in this 
paper, we will present an explanation to these questions from another point of view, 
namely using a dual 2 + 1-dimensional Maxwell electrodynamic (MED) description, 
and we will show that the main features observed can be understood naturally in this 
framework, although qualitatively in some aspects. 

This dual electromagnetic description of the interacting bosons model was 
obtained by the duality transformation early in the study of Helium superfluidity, [H] [H] 
and then in the past decades it has been discussed in the contents of superconductivity 
and superfluid film by various authors. [TT)||11| Compared to the model discussed 
before, we would like to emphasize some important differences of the system studied 
here. First, a large amount of angular momentum have been imparted into the atoms 
and the rotational invariance is restored before the formation of condensate, therefore 
the total angular momentum of the condensate should be conserved to a non-zero 
value. The second, the atom removal laser violates the current conservation condition. 
And the third, the system is dilute and the interaction between atoms is weak, the 
sound velocity is consequently relatively small. We will show that these differences 
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lead to these exotic phenomena in these experiments, in other words, these observed 
phenomena reflect some intrinsic properties of the rotating BEC. 

2. The Dual Picture 

In this section we will firstly make a brief description of the dual picture for the 
interacting bosons system, for the details we refer to the Ref.[§][TU]|ll|. Beginning with 
the coherent state path integral formulism. |12j we write down the quantum mechanical 
propagator 



Z = J D[il>,ip]exp{i I dtd 2 fC}, (1) 
where the Lagrangian 
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g is the interaction strength and equals to ATrh 2 a sc /m for ultracold atomic system, 
where a sc is the s-wave scattering length between atoms. Here ip is a complex field 
describing the bosonic many-body system, and ip is its complex conjecture. We can 
factor ip as y/pe l6 V 9 " . The field p describes the density of bosons. The phase of ip 
contains two parts, 9 is the smooth part free from singularity and 6 V is that caused 
by the presence of vortices. 

First of all, we can introduce an auxiliary field J by performing a Hubbard- 
Stratanovich (H-S) transformation, the physical meaning of which is the superfluid 
current. By defining the three vectors J = (hp,hJ) and V = (9t,V), in the low 
frequency and long wavelength limits one can neglect the second order gradient term 
\Vp\ 2 / p, and obtain that 

Z = J D[J]D[p]D[0]D[0 v ] exp U J dtd 2 f (-J • V0 - J ■ V6 V + ^|J| 2 - |p 2 

Integrating over 9, one can naturally obtain a 5-function constraint V • J — 0, which 
is just the conventional current conservation condition 

| + V.J = 0. (4) 

Owing to this constraint, we can introduce a three-component gauge field b defined 
as (po,b), which satisfies 

J = V x b, (5) 

the density and the current of the bosons are thereby determined by the gauge field. 
Hence, p = V x b is compared to a magnetic field perpendicular to the 2-dimensional 
plane. We can expand p as p + 6p(r, t), p is the average density which is compared to 
an external magnetic field, while 8p(f, t) is the density fluctuation which is compared 
to the electromagnetic wave. The electric field E is defined as V6o — dob, and it is 
perpendicular to the boson current J. Therefore the Eq.Q can be rewritten as 

— + Vx£ = 0, (6) 
which is just the Maxwell equation in a monopole free space-time. 



(3) 
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With replacing the p by p in the denominator, the last two terms in Eq.JJJl can 
also be reexpressed in terms of the gauge field 6, which is just the dynamic term 
_pnvp f t ne gauge field as in the conventional MED. ft determines the dynamic 
properties of boson density fluctuation, i.e. phonon, and recovers the well-known result 
for the sound velocity s — \J gp/m, which plays the same role as the light velocity in 
2 + 1-dimensional MED. 

The coupling between the vortices and bosons can be illustrated as 

-J-VB v = -b-(vx.V6^\=-b-J v . (7) 

Here J v is defined as V x V0„ . It can be verified that the zero-component of J v , which 
is denoted by p v , represents the vortex density, and the other two spatial components 
represent the vortex current. The Eq.Q indicates that the vortices can be viewed as 
charged particles, which are coupled to an electromagnetic field via gauge coupling. 
The coupling constant, that is the charge q of the charged particles, is unit. In the 
Coulomb gauge, the scalar potential 6o is instantaneous, one can integrate it out 
and obtain the mass of vortex m„, which is equal to ph 2 In ^/(47rs 2 m) for a uniform 
condensate, and the mutual logarithmic interaction between vortices described by 

dnt = J d 2 rd 2 Pp v (r) In ^^p„ {?). (8) 

Here R is the radius of the condensate and A is the radius of vortex core. 

Thus we have obtained a dual electrodynamic description of quantized vortices, 
which is briefly summarized in the following table. In the dual picture the vortex 
lattice state corresponds to the Wigner crystal state of 2-dimensional electrons in the 
external magnetic field. 



Original Picture 


Dual Picture 


vortex 


charged particle 


average 
superfluid density p 


external magnetic field 


superfluid current J 


electric field = e^ v J v 


density fluctuation 


electromagnetic wave 


sound velocity 


light velocity 


the coupling between vortex 
and sound wave 


gauge coupling 



Such a description can be directly applied to the confined BEC when it is in 
the Thomas-Fermi regime. The average density here is therefore determined by the 
trapping potential, as well as the centrifugal repulsive force caused by rotating, and we 
will take it as inverted parabola in this regime. Here we would like to point out that 
the main assumption leading to the electrodynamic description is that the vortices 
are treated as point-like particles, the physics inside the vortex core is not considered, 
as a consequence we take the mass of vortex as its electromagnetic mass throughout 
this paper [TT]. Besides, in the following treatment we view the system as quasi-two 
dimensional and neglect the density inhomogeneity when estimating the vortex mass. 
In the quantum Hall mean field regime[T5|. it is shown that the physics of vortex core 
will become important. However, before entering this regime, and when the Thomas- 
Fermi theory is still valid, these neglected effects will only slightly renormalize the 
vortex mass as well as the interaction between vortices. 
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3. Explanation of the Experiments 

3.1. Vortex Number Increasing 

In JILA experiments the rotational invariance is restored before condensation, 
therefore the condensate should obey the total angular momentum conservation 
constraint. Such a system can be described by adding a Lagrange multiplier term 
to the Lagrangian. 

Z = J d^pdnexp{i J d 2 rdt£(?P,iP) - n($LiP - L a )}. (9) 

It is obvious that one will obtain a (5-function which means that the total angular 
momentum must always equal to Lq after integrating f2 out. The physical meaning of 
Q is in fact the effective rotation frequency. Following the procedure discussed in above 
section, it is found that only three additional terms, namely — 2mbon+^mil 2 pr 2 +flLo, 
will be added to the Lagrangian. One can see that only the zeroth component of the 
gauge field responds to the rotation. Integrating f£ out and neglecting the density 
fluctuation, we obtain 

(-2mb + L ) 2 

C = — ; — h other terms. (10) 

2m (pr 2 ) 

Following the same step which leads to Eq.©, integrating out bo, we can obtain a 
local chemical potential term in the following Lagrangian 

Ant = - / d 2 rd 2 ?p v (r)G(\r-?\)p v {?) + / d 2 r^ ( [ d 2 r'G(\r ~ r'\) ) p v (r) 



(pr 2 

£ mt + I d 2 rp(r)p v (r). (11) 



Here G(\r* — r"|) is the Green's function for the &o field acquiring a very small mass 
2m /(pr 2 ), which will vanish in the large atom number limit. The first term in 
the Eq. I|ll|l indicates that the vortex system can be viewed as an interacting 
Coulomb gas system, and the second term defines a local chemical potential p(r) 
for vortices emerging from the angular momentum constraint, which is proportional 
to Lo/(pr 2 ). The atom removal laser, which removes atoms with lower angular 
momentum, decreases the atom density p with little change to the total angular 
momentum Lq, leads to the increase of the vortex chemical potential. Notice that 
the interaction energy £i nt depends quadratically dependence on vortex density p v , 
and the interaction between vortices are repulsive, so 

da 8 2 £i n t , . 

op v 5 pf, 
This shows the increase of vortex density. 



3.2. Aggregate 

For simplicity we assume that the atom removal laser acts only at the center point. 
Extension to the Gaussian shape laser is straightforward and the qualitative results 
will not be changed. When the laser is applied, a centripetal atom current will be 
induced and the current conservation condition will be violated at the center. In the 
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Figure 1. (a) Atom removal laser induces a current of atom in the original picture 
(b) Vortices are viewed as charges moving under a perpendicular magnetic field 
in the dual picture. The thin lines represent the magnetic and electric field, the 
thick lines represent the atom removal laser or monopole current. 



dual picture, a circular electric field will be generated and the Maxwell equation should 
be modified in the following way 

^L- + WxE = I5(r). (13) 

Here the value of I equals to the number of atoms removed from the condensate per 
unit time. Therefore it is modelled as a monopole current running through the plane. 
The magnitude of the electric field is approximately J/(27r|r|) nearby the center. For 
the case of continuous laser the vortices as charges will be accelerated azimuthally by 
the electric field, and the Lorentz force will drive the vortices to aggregate toward the 
center, and eventually to form a giant vortex. 

3.3. Two Kinds Oscillation 

After a pulse laser which only acts for a short time AT, the vortices will acquire a 
velocity v and leave the balance positions of the lattice. For the vortices closest to the 
center, the magnitudes of their velocities can be approximated by 

v = EqAT/m v ~ IATh/(2Trm v d), (14) 

d is the vortex lattice spacing. Therefore, there are two forces which can cause the 
oscillation motion, one is the Lorentz force /i , and the other is the mutual interaction 
force between vortices f^. It is easy to obtain the ratio of the two forces: 

ll = /AT—. (15) 
J2 m v 

For these experiments we can naturally assume that the number of atoms removed 
from the condensate / is an increase function of the power of removal laser. Thus 
Eq. l|15|) indicates that the Lorentz force is dominant for a strong laser pulse, such as 
the pulse used in Ref.p^. In this limit, the vortices will execute circunmutation. When 
the vortices move very close to each other, they will merge into a giant vortex, and 
separate into several single vortices again when they oscillate away from each other. 
The mutual interaction force will be dominant for the weak pulse case. In this limit, 
it will cause a collective oscillation of vortex lattice known as Tkachenko oscillation. 
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The oscillation frequency u> depends on the initial conditions of the condensate, 
such as the number of atoms and the rotation frequency f2. In the strong pulse limit, 
for the motion is mostly determined by the Lorentz force, w = Tip/m v . Because the 
density p is proportional to the number of atoms N and is a decreasing function of 
the rotation frequency f2, the oscillation frequency is linearly dependent on N and 
increases with decreasing. 

In the weak pulse limit, the oscillation frequency u> is proportional to yfp/d. It is 
not obvious to tell whether ui will decrease or not when the condensate rotates faster, 
because both p and d will decrease with the rotation frequency f2 increasing. We make 
a rough estimate as following: For a fast rotating BEC confined in a harmonic trap 
with trapping frequency u>±, the effective confining potential is {uj\ — Q 2 )r 2 /2. Using 
the Thomas-Fermi approximation we can obtain that the condensate density nearby 
the center p is proportional to ^Jl — Of, and then using the formula d 2 oc l/p v — 
h/(2m£l), we can obtain that 



C 



\ 




(16) 



Here C is a complicated constant, which is related to the atom mass, vortex mass, trap 
frequency and the density distribution along the third spatial dimension. The Eq. i|ltjfl 
shows that ui is a decreasing function of f2 when Q / co± > 1 / V2, which is shown in the 
left side of Fig.©[Hj. 




0.86 0.88 0.9 0.92 0.94 0.96 0.98 

Omega(in the unit of trap frequency) 



Figure 2. The oscillation frequency of vortices u as a function of rotational 
frequency of condensate f2(in the unit of uj±). 



From above discussion we can see the atomic current toward the center is of 
essential importance for causing oscillations. Hence if there is any other ways to 
generate a similar atomic current, such as using optical dipole force to draw atoms 
into the middle of the condensate in the Ref.[2], the similar oscillation of vortex lattice 
was observed. 



3.4- Damping Effect 

One feature of these oscillation motions shown in these experiments is that the 
oscillation amplitudes are heavily damped, especially for the Tkachenko oscillation 
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Figure 3. Heavily damped oscillation amplitude as a function of time t(in the 
unit of oscillation period T) 



reported in Ref.|2]- It is quite natural in the view of the dual picture. Because 
the coupling between the vortex current and the phonon field is a gauge coupling, 
an accelerated vortex will radiate sound wave in the same way as an accelerated 
electron radiating electromagnetic wave. In the view of hydrodynamic picture, when 
a vortex is accelerating, the circular current J, as well as the singularities of field J, 
can not be stationary in any inertial frame. Due to the current conservation condition 
Eq.J3J, a time dependent density fluctuation Sp(f,t) will be generated. Hence the 
oscillating vortices will induce sound wave, lose their energy and result in the damping 
of amplitude. In the dual picture this process can be calculated through a standard 
MED calculation. 

The retarded Green's function for 2 + ID MED is 

And the gauge potential b ra d for the radiation field is 

t-\f-f(t')\/s 

*> , V{t ' ] (18) 

y/s 2 (t-t') 2 - (f-r(i')) 2 

Here V(t') is the vortex 3-current. Unlike the 3+1D case, the Green's function contains 
the step function instead of Dirac's 6 function, so one can not easily find a general 
expression for the radiation energy of a charge with arbitrary acceleration. |15 j |16 | Here 
we only consider two special cases, the linear oscillation V(t') = (v cos(wt'), 0), and 
circular motion V(t') = (— v sin(wi'), v cos(urf')). Following the approximation made 
in Ref.^nji we replace \r — r(t')\ with R and calculate the energy radiation along large 
circle \r\ = R. Then we can perform an analytical calculation for the average energy 
radiated over one period, 

T 2tt 

P = i J J dtdOra-r, (19) 
o o 



where a is the Poynting vector. For these two special cases, the result is 

P =f4 = i?pL. (20, 

2ms z zg 

The radiation rate is inverse proportional to the interaction strength g, so the damping 
effect is remarkable for the weak interacting boson gases. Recalling the approximate 
expression for the vortex mass m v = ph 2 ln(i?/A)/(47rs 2 m), Eq. 120fl can be reexpressed 
as 

1 pk 2 2 1 47TW o 

± = — — — ojv — — wm v v = ttE. 

2s 2 m 21nf In f 



Hence 

dE Attlu 



dt ~ In | 



E, (21) 



this indicates that the oscillation energy exponentially decays as exp{— t/r}, and 
the ratio of the exponential damping time r to the oscillation period T is equal 
to ln^/(87r 2 ). A typical oscillation curve in the current experiment conditions 
is schematically shown in the right side of FigEI one can see that the oscillation 
amplitude will decay rapidly within one period, and this is qualitatively coincident 
with the experimental observation However, there are also a little quantitatively 
disagreement between the theoretical damping curve and the experiments, in Ref.pQ 
the amplitude vanishes within about three periods and in Ref. jzj it is about one and a 
half periods. This disagreement occurs because in the real case the vortex mass should 
be modified due to the spatial density inhomogeneity, the third spatial dimension effect 
and the edge effect. The mass of vortex, which is proportional to the self-energy, is 
therefore changed with the position of vortex. ^7] It is expected that the damping 
curve will agree better with the experiments when these subtle effects are carefully 
considered. 



4. Summary 

According to the effective Maxwell electrodynamic model the vortices in the BEC can 
be viewed as an ensemble of classical charges moving in an external magnetic field. So 
far we have applied this description to understand the observed phenomena of vortices 
induced by selective removal atoms from a fast rotating BEC. The atom removal laser 
applied to the center of the condensate will cause two leading effects in the view of dual 
picture. One is the increase of the vortex chemical potential and then subsequently 
the increase of vortex density. The other is that it will induce a circular electric field, 
and each vortex will acquire an azimuthal velocity. We clarified that there are two 
forces that can cause oscillation motion, and discussed the the relation between the 
oscillation frequencies and the number of atoms together with the effective rotation 
frequency. We showed that the dispersion mechanism for the vortices oscillation comes 
from the emission of phonon, which is the same as the radiation of the electromagnetic 
wave of an accelerated charge. And the damping effect is more remarkable in these 
experiments because the interaction strength is weak. 

In this paper the dual description of vortices has been used to explain current 
experiments. Although it is not yet a detailed quantitative description, it provides a 
natural physics picture of most key experimental observations, which convinces us to 
believe that this description is valid to the extent that the BEC is fast rotated but 
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still not entering the quantum Hall mean field regime. We hope this will start a new 
route to study the issues in the rapidly rotating BEC system. 
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